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Abstract 

We reexamine the strength of the first order phase transition in the electroweak theory 
supplemented by an extra Higgs doublet. The finite-temperature effective potential, V e s, is 
computed to one-loop order, including the summation of ring diagrams, to study the ratio 
4> C /T C of the Higgs field VEV to the critical temperature. We make a number of improve- 
ments over previous treatments, including a consistent treatment of Goldstone bosons in 
V e f[, an accurate analytic approximation to V e R valid for any mass-to-temperature ratios, 
and use of the experimentally measured top quark mass. For two-Higgs doublet models, we 
identify a significant region of parameter space where 4> C /T C is large enough for electroweak 
baryogenesis, and we argue that this identification should persist even at higher orders in 
perturbation theory. In the case of the minimal supersymmetric standard model, our results 
indicate that the extra Higgs bosons have little effect on the strength of the phase transition. 



1 Present address: Department of Physics, UCLA, 405 Hilgard Ave., Los Angeles, California 90095-1547 



1 Introduction 



There is now convincing evidence that the baryon asymmetry of the universe could have been 
created at the electroweak phase transition (EWPT) [1], with a minimal amount of new, yet 
plausible physics beyond the standard model [2]- [4]. However for electroweak baryogenesis to 
work, it is necessary that the baryon-violating interactions induced by electroweak sphalerons 
be sufficiently slow immediately after the phase transition to avoid the destruction of the 
baryons that have just been created. This condition is fulfilled if the vacuum expectation 
value (VEV) of the Higgs field in the broken phase is large enough compared to the critical 
temperature at the time of the transition (see for example ref. [5] and references therein), 



In the standard model, the bound (1) holds only for very light Higgs bosons [6]- [10] already 
ruled out by experiment, m^o < 20 GeV, according to a recent nonperturbative study of the 
phase transition [11]. In addition to the difficulties with producing a large enough initial 
baryon asymmetry, the impossibility of satisfying the sphaleron constraint (1) in the standard 
model provides an incentive for seeing whether the situation improves in extended theories. 

The EWPT in two-Higgs doublet models has been previously studied in references [12]- 
[14] . In the present work we improve on the approximations made in the previous studies in 
the following ways. 2 

• Avoidance of the unitary gauge, which is known to be problematic: it has been shown that 
it is necessary to resum an infinite number of contributions in the perturbative expansion 
to make the determination of finite-temperature quantities agree between unitary gauge and 
renormalizable gauges [15]. 

• Inclusion of the /i 2 .*^^ term in the Higgs potential necessary for CP violation, which is 
crucial for the baryon asymmetry production mechanism. In contrast to ref. [14], we find 
that the strength of the transition does not necessarily grow with increasing /x 2 .. 

• Use of the ring summation [6]- [8], [16] for all the bosonic thermal loop contributions to the 
effective potential. Furthermore we compare two different prescriptions for implementing 

2 Thc list is not meant to imply that all of the previous studies were subject to all of the limitations 
mentioned here. For example, only ref. [13] used unitary gauge. 



C /T C > 1. 
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ring-improvement. The differences are indicative of the theoretical uncertainties inherent in 
using one-loop perturbation theory. 

• Use of a critical temperature T c more closely corresponding to the onset of bubble nucle- 
ation, namely the temperature where the two minima of the potential are degenerate, rather 
than that where the second derivative of the potential vanishes at the origin of field space. 
We find a quantitative difference between the two methods when determining whether a 
given set of parameters in the model is really consistent with the sphaleron bound. 

• Use of the experimentally measured top quark mass, 175 ± 6 GeV [17], which was not 
available when the previous papers [12]- [14] were written. The large mass of the top quark 
greatly reduces the strength of the phase transition in the standard model [8]-[10]. 

• Estimating the expansion parameter which quantifies the convergence of the perturbative 
expansion at the phase transition, near the critical value of the Higgs field VEV[18, 8]. This 
is the major potential weakness of the perturbative approach being used here. However we 
find that according to this criterion, most of the promising regions of parameter space are 
relatively safe from the infrared divergences that would cast doubt on the convergence of the 
loop expansion. 

• Inclusion of the Goldstone boson contributions to the effective potential. These can become 
important when the lightest Higgs boson h° starts to become as heavy as the other particles 
in the theory. Although the case of large m h o is uninteresting in the standard model, from 
the point of view of having a strongly first order phase transition, in extended models the 
effect of large m^o can be compensated by new particles. 

• Avoiding the expansion in mass over temperature that would limit the scope of our inves- 
tigation. This turns out to be quite important since the region of parameter space where 
the phase transition is strongly first order is largely where the conventional expansion of the 
effective potential in powers of M/T is starting to fail badly. 

A price we pay for these improvements is the restriction of the potential such that the 
ratio of the two Higgs VEV's (tan (3) is unity, to make the problem more tractable; thus we 
are not exploring the entire parameter space of the class of models of interest. On the other 
hand, this reduction has the advantage of making it easier to characterize the favored values 
of the parameters in the Higgs potential. Furthermore a number of studies of the minimal 
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super symmetric standard model (MSSM) indicate that the strength of the phase transition 
is generally greatest for small values of tan (3 [19]- [24]. In section 4 we discuss the significance 
of our results for the MSSM in the region of tan/3 = I. 



2 Effective Potential 

2.1 Zero temperature; Goldstone boson contributions 

We begin by constructing the effective potential for the CP-even component of the Higgs field 
that is responsible for electroweak symmetry breaking. At one loop and at zero temperature, 
each particle makes a Coleman- Weinberg- type contribution to Kff,o(0) [25]. It can be written 
in the form of a supertrace 3 involving the matrix M of all the field-dependent particle masses, 

W0) = \{<f ~ v 2 ) 2 + + ^Str M 4 (0) (log _ , (2) 

where the pieces proportional to A and log(yU 2 ) are counterterms that must be fixed by 
renormalization conditions. A convenient choice of the latter is that the tree-level definitions 
of the VEV (v = 246 GeV) and the Higgs mass m 2 = V"(v) = 2\v 2 should be maintained. 
These conditions are given by 

1 / M 2 \ 

=V'(v) =Av + —StrM 2 M 2 'l\og—-lj (3) 

(M 2 M 2 " + (M 2 ') 2 ) (log % - l) + (M 2 ) 



<o =V"(v) = 2Xv 2 + A + T^Str 



(4) 



where the prime denotes d/d(j) and all masses are evaluated at (f) — v. One can always choose 
A such that eq. (4) is satisfied when m 2 h0 = 2\v 2 , so the tree-level relation between m^o and 
the VEV is preserved. 

However there is a technical problem with eq. (4) : in the Landau gauge which is the most 
convenient one because of the decoupling of the ghosts, the Goldstone bosons are massless at 
= v , but they have a nonvanishing value of M 2 '(v), giving a logarithmic infrared divergence 
due to the term (M 2 ') 2 logM 2 , which does not occur in other gauges. Yet physical masses 
must be gauge-invariant and IR-finite. The problem is that m 2 h0 as defined in (4) is not the 



3 each real scalar field has a weight of +1 in the supertrace, and each Dirac fcrmion has a weight of —4 
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physical pole mass, which must be evaluated at an external momentum of p 2 = m 2 , but 
rather the off-shell self-energy evaluated at p 2 = 0. To avoid the IR divergence and the gauge 
dependence we should really evaluate the pole mass. By computing the Feynman diagrams 
of the type shown in figure 1 which are responsible for the IR problem, it is easy to see how 
to fix it: in eq. (4) one should make the replacement 



log M 2 -> Re J dx log (M 2 - x(l - x)m 2 h ^j 



log R + 2V4i? - 1 tan- 1 ( 7 =f) , R > 
= F(M 2 ,m 2 ); R = M 2 /m 2 1 o, 



i. 

4' 



Using the prescription of eqs. (4-5), and the convention m 2 = 2Xv 2 , it is straightforward to 
solve for the renormalization constants: 

log /i 2 = Str [(m 2 M 2 '/v + (M 2 ') 2 ) F(M 2 , m 2 ) - (M 2 M 2 '/v) log(M 2 )] /Str (M 2 ') 2 ; 

We used the fact that for any renormalizable theory, M 2 " = M 2 ' /v, to simplify the result. 

By this procedure we are able to consistently include the contributions of the Goldstone 
bosons into the effective potential (2). In the standard model, for small values of m 2 h0 which 
were of interest for maximizing the strength of the phase transition, both the usual Higgs 
boson and thus the Goldstone bosons make numerically small contributions to the effective 
potential and can thus be safely ignored. In two Higgs doublet models however, we will show 
that it is possible to get strong phase transitions even when m^, becomes rather large, and 
therefore the Goldstone boson contributions can be nonnegligible. 



O 



Figure 1: self-energy diagram responsible for IR-divergent contribution of Goldstone bosons 
to the Higgs boson mass at p 2 = 0, in Landau gauge. 
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2.2 Finite temperature 

At one-loop order but now at finite temperature, there is an additional contribution to the 
effective potential given by [26] 

7^4 poo / \ 

AV eS , T (cf>) = ^jStrjf dxx 2 \n (l T e -V* 2 +^W 2 j (7) 

where the sign is — for bosons and + for fermions. For small M/T the contribution per 
degree of freedom can be expanded as 4 

T/ _ w 2 T 4 | M 2 T 2 M 3 T M 4 (] na .(M 2 \ _\ 

, M¥W-MM' (a-3)!!((2i-i) bosons- 

+ 2 Z^V 471 " 2712 / (2/)!!(Z+l) ' UUSUI1S, 
T/ /„\ _ 7tt 2 T 4 , M 2 T 2 , M 4 f-i (M 2 \ n \ 

V s ,i[n) - — + + (log [yrj-Cfj 

n , 

m 2 t 2 v- f —M 2 \ (2i-3)iic(2f-i) ( 2 2i-i A fermions- 
2 2^\4n 2 T 2 ) (2/)!!(/+i) { z L ) ' iermions, 

i=2 

c 6 = 3/2 + 2 log 4tt - 2 7b 5.40762; C/ = c b - 2 log 4 2.63503 (8) 

where ( is the Riemann £ function. The term cubic in the boson masses is especially impor- 
tant because it gives rise to the barrier in the potential that makes the transition first order, 
as needed for baryogenesis. For large M/T on the other hand, the contribution from either 
bosons or fermions has the asymptotic expansion 5 

For the purpose of numerical analysis of the phase transition, it is useful to have an 
analytic approximation to the exact expression (7) since it is computationally expensive to 
evaluate the integral. However, by smoothly joining the small M/T expansion (8) with that 
for large M/T, eq. (9), one can obtain an excellent approximation to the exact integral. 
To do this optimally, one should choose the value of M/T where the derivatives of the 
two expansions match each other as the transition point. Using VJ(3) for the large M/T 
approximation, we found good matching to the exact integral by using V^b(3) for bosons 

4 see eq. (3.32) of ref. [8] 
5 see eq. (B5) of rcf. [9] 
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and V^f(5) for fermions. This gives an approximation with a relative error which is less than 
0.5% for M/T = oo, and negligible for small M/T. Specifically, we took 



A\/ cff ,T = Tr 



[Q(x b - (M/T) 2 ) V B)b (3) + 0((M/T) 2 - x b ) (V,(3) - 6 b )\ P b 

[Q(x f - (M/T) 2 ) K,f(4) + 6((M/T) 2 - x f ) (VJ(3) - 5/)] P, 1 , (10) 



where 6 is the step function with x b = 9.47134 and Xf — 5.47281, and P b j denotes the 
projection operator for bosons or fermions, respectively. The small constant shifts of Vj(3) 
are made so that the function as well as its derivatives match at the transition point: 5 b = 
3.1931 x 10~ 4 and 5 f = 4.60156 x 10" 4 . 

The above discussion implicitly assumed real- valued masses, but for completeness there is 
one further matching of large and small M/T behavior which must be considered separately: 
it can happen that some of the bosons have large and negative values of M 2 , corresponding 
to imaginary M/T. While it may be arguable whether eq. (7) has a meaningful physical 
interpretation for negative M 2 , to the extent that it has any meaning, it is clearly necessary 
to use the large M/T expansion when \M 2 /T 2 \ is large. For negative values of (M/T) 2 
the matching conditions are different; in eq. (10) one should use x b = —6.84368 and 5 b = 
0.045374. Since large values of |M 2 /T 2 | always imply the decoupling of the heavy particles, 
none of our results concerning the interesting regions of parameter space, where the sphaleron 
constraint eq. (1) is satisfied, will depend on the precise value of the effective potential when 
the squared masses are large and negative. 

One can easily improve the above one-loop result by resumming a subclass of thermal 
loops known as the ring diagrams, which amounts to finding the thermal corrections to 
the boson masses 6 due to AV eStT {<P), typically of the form M 2 (0, T) = M 2 (0) + g 2 T 2 , and 
substituting them back into the expression for the total effective potential [16]. In fact the 
choice of exactly how to resum is not unique, and Arnold and Espinosa [8] have justified 
the consistency of the widely-used procedure of replacing M 2 (0) — > M 2 (0, T) only in the 



6 the reason this does not apply to the fermion masses can be understood most easily in the imaginary 
time formalism of finite-temperature field theory, where the effective squared masses of the Matsubara modes 
are M 2 {<j), T) + (2irnT) 2 for bosons and M 2 (0, T) + (27r(n + ±T) 2 for fermions. Only for the n = modes of 
the bosons can there be an infrared divergence due to vanishing M 2 (<j)) which would make it important to 
include the perturbative g 2 T 2 contribution to M 2 . 
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cubic term M 3 T of eq. (8) rather than everywhere in the effective potential. Defining the 
unresummed total one-loop effective potential as 

V tab = V eS , + AV eB ,T, (11) 
the ring-improved potential in the two methods can be written as 

King = V tot [M 2 ((f), T)] , Parwani method 
King = V tot [M 2 ((j))] + — ( m3 (0) - M 3 (0,T)), Arnold-Espinosa method 

bosons 

We shall investigate and compare both of these prescriptions. Since they differ by terms 
which are of two-loop order, they can give us some idea of the uncertainties in our calculation 
coming from the neglect of higher orders in perturbation theory. 

The ring summation is essential for correctly estimating the contributions from the lon- 
gitudinal gauge bosons, which get a much larger thermal mass than their transverse counter- 
parts, because this tends to reduce the effectiveness of the cubic term in making the phase 
transition first order. In addition, the thermal corrections decrease the region of field space 
near = where some of the squared masses of the bosons become negative, giving rise to 
a complex potential. However the imaginary part of the potential can still persist because 
of the genuine instability of field configurations with small values of 0. We find that the 
real part of the potential often exhibits an additional local minimum especially when some 
of the bosons have a small or negative M 2 , as illustrated in figure 2. Such a situation indi- 
cates that the negative squared masses of the Higgs bosons are playing a large role and that 
perturbation theory should not be trusted due to infrared divergences in the effective three 
dimensional theory which describes the high temperature limit. We discuss the consequences 
of this phenomenon further below. 

3 Two Higgs Doublet Model 

Let us now introduce the new physics. For the two Higgs doublet model, including the 
couplings to the top quark, the potential is taken to be 

K-doublet(<f>l, $2) = - £®\®2 " Vl®\®2 ~ l4'&2*l + y I)' + y ( $ 2$ 2? 
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Figure 2: Evolution of the effective potential with temperature, showing the pathological 
case of two nontrivial minima, indicative of negative squared boson masses. The parameter 
values are m hQ = 125 GeV, ji\ = (30 GeV) 2 and all other Higgs boson masses 187 GeV. 

+ / il ($ t 1 $ 2 )(^$ 1 ) + /i 2 ($ t 1 $ 1 )($^ 2 ) + h 3 (($I$ 2 ) 2 + ($^i) 2 )) 

+ T L (y 1 $ 1 + y 2 $ 2 )t R + c.c. (12) 

It incorporates a softly-broken discrete symmetry, $j — > — $j, in order to reduce the number 
of possible couplings. This is the same model as discussed in ref. [13], except for the addition 
of the term /i|${$ 2 , which is necessary for having CP violation in the Higgs sector. To 
simplify the analysis, we follow ref. [13] by setting fif = //|, Ai = A 2 and y\ = y 2 - Furthermore 
we assume that //| is real and positive. 7 This choice insures that the symmetry will break 
along the direction $i = $ 2 = |(0, <p) T at any temperature. The resulting potential for <p is 
identical to the standard form A(0 2 — t> 2 ) 2 /4 used in (2) if one identifies 

A = i(A 1 + /i 1 + /i 2 + 2/ i3 ); v 2 = d±A, (13) 

The spectrum consists of the light Higgs boson h° associated with the field 0, three Goldstone 

bosons G°,G ± , the charged particles H ± , and the neutral scalar and pseudoscalar, H° and 

7 The sign of n\ is conventional since it can always be changed by the field redefinition $i — > — $i. 
Moreover from the results of ref. [14] it appears that the strength of the transition depends only upon \fi^\ 
when ^3 is complex. 
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A repsectively. 

All that remains to complete the definition of the effective potential is to specify the field- 
and temperature-dependent masses. For the gauge bosons these are the eigenvalues of the 
mass matrix in the basis W ± , W 3 , B, 

\ / 9 2 



m 2 (0,T) = 



V 



9' 



9* 
99' 



99 

g' 2 



J 



g' 



V 



\ 



g' 2 ) 



2, longitudinal; 

0, transverse. 



(14) 



Notice that only the longitudinal gauge bosons get a thermal mass at this order (for the 
transverse bosons there is a magnetic mass of order g 2 T which we neglect) and the coefficient 
2 in the thermal correction would be 11/6 in the absence of the second Higgs doublet. Also 
the longitudinal photon is no longer massless at finite temperature, due to mixing between 
B and W 3 . The top quark mass is 





mt(<f>,T) = m t - 



(15) 



Although the dispersion relation for the top quark does get temperature corrections which 
can be interpreted as a thermal contribution to its mass squared, these do not appear in the 
resummation of the ring diagrams because it is a fermion, as explained above. For the Higgs 
bosons, it is convenient to eliminate all the dimensionless couplings and express the general 
masses in terms of the zero-temperature ones and the parameter n\ [4]: 




-m 



h° 



(m 



-m 



h° 



4fi 



h° 

2 4 

3 ^ 



2m 2 ) € 



V1V2 \ rp2 

1 



[m 



h° 



h°,G°,G ± ; 
A , H°, H ± ; 



A , Sg 2 + g' 2 yj + yj 
a = 2 + ^^ + ^ 



-4^ + 2mf)$) + (a-^)T 2 , , 
j2 { 4A + ^(-^3 + 2m ^± + m 2 H0 + mio)} ; 



(16) 



For comparison, the standard model can be recovered by ignoring the masses of A°,H°, H* 1 
and the term in curly brackets in a, and setting y 2 = 0. The particle multiplicities are 1 
for each of the Higgs bosons and the longitudinal gauge bosons, 2 for each transverse gauge 
boson, and 12 for the top quark. 
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4 Phase Transition Results 



In studying the quantity <p c /T c one should in principal take care to distinguish between 
various definitions of the critical temperature. The T c relevant for baryogenesis is the tem- 
perature when tunneling from the false vacuum to the true vacuum takes place, after which 
the phase transition quickly completes. To determine it precisely, one demands that the 
tunneling probability becomes sufficiently large compared to the Hubble expansion rate, 
which requires solving for the bubble configurations and computing their action [10]. In the 
standard model, rather close upper and lower bounds on T c can be obtained because it is 
always above the temperature To where V"(0,Tq) = 0, and below the temperature T\ when 
the nontrivial minimum v(T) in the potential becomes degenerate with the minimum at the 
origin, and these temperatures are all within a few tenths of a percent of each other. However 
in extended models it might be important to make the distinction between the various crit- 
ical temperatures because they can differ more substantially, and also because v(T)/T can 
change rather quickly near the critical temperature. In the following we compare the results 
of using either definition of the critical temperature and find some noticeable differences, 
particularly for small values of the Higgs mass parameter //|. 

Because the masses of the heavy Higgs bosons (A , H° and H^) have such a similar 
form in eq. (16), we have simplified our scan of the parameter space by taking them all 
to be equal to each other. Figure 3 shows the contours of constant <f)/T in the plane of 
m A o = m H o = m H ± versus m h o, using three different methods. The first column uses 
T as the definition of the critical temperature, and the Arnold-Espinosa method of ring- 
improvement. The second column is the same except using Ti for the critical temperature. 
The third column also uses 7\, but the Parwani ring-improvement prescription. Each row is 
for a separate value of the potential parameter //g, ranging between and 6 x 10 4 GeV 2 . 

The regions in white are those with C /T C < 1, and therefore indicate the parameters 
for which electroweak baryogenesis would not work in this model. In all three methods 
one sees that for a given value of m^o, there exists a range of the other masses such that 
4> C /T C > 1. When the other masses are greater than this range, it turns out that the 
nontrivial minimum is higher than the one at = even at zero temperature, meaning that 



10 



To 



Arnold-Espinosa 



Parwani 



m A o,HO,H^ 6 



t4 



(4 = 6 




Figure 3: Contours of constant <fi c /T c (first three columns) and the perturbative expansion 
parameter e (defined in eq. (19)) for three values of in the plane of m^o = m#o = m#± 
versus m h o. The first three columns use different approximations for the critical temperature 
or the effective potential, described in the text. Units are 100 GeV for masses, (100 GeV) 2 for 
/i§. Increasing values are represented by darker shades, with white being < 1 and increasing 
in steps of 1. The regions labeled by "m." in the last column are where the false vacuum 
is metastable, and those labeled "s.o." are where the phase transition is second order, as 
illustrated in figure 4. In neither case is there a first order transition. The values of e are 
computed using the boson masses corresponding to the second column. 
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Figure 4: Examples of the evolution of the potential with temperature which illustrate the 
absence of a first order phase transition. In the first one the false vacuum is metastable, and 
in the second the phase transition is second order. The parameter values are = 60 GeV, 
other masses = 325 GeV, ji\ = for the first one, and m h o = 300 GeV, other masses = 250 
GeV, ^1 = for the second. 

the electroweak symmetry breaking vacuum is metastable and there is no phase transition, 
unless supercooling into the metastable minimum occurs. On the other hand, when they lie 
below the allowed bands, the transition becomes second order, because there is never any 
barrier separating the minimum at <fi = from the symmetry-breaking minimum. In either 
case we define <fi c /T c to be zero. Such examples are illustrated in figure 4. 

However the three methods of constructing the finite-temperature potential differ as to 
the precise values of C /T C when /x| is small, with the Parwani ring- improvement giving larger 
values than that of Arnold and Espinosa. For large /i 2 , all three give similar results. Where 
the two methods of resummation differ is indicative of where perturbation theory may be 
unreliable. There is another curious phenomenon for small but nonzero values of for 
example ji\ = (30 GeV) 2 , which is not shown in figure 3: the region allowed by baryogenesis 
appears to be significantly enlarged in the first two methods (To and Arnold-Espinosa), 
though not in the third. Examination of the potential shows that the region of enlargement 
is always accompanied by the existence of double minima, as in figure 2. Therefore the 
augmentation of the baryogenesis-allowed parameter space in these cases is spurious. This 
is one possible explanation for our failure to corroborate the claim of ref. [14], which says 
that increasing /i| generally has a favorable effect on the strength of the phase transition. 

In addition to comparing different methods of ring-improvement, we can also quantify 
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the convergence of the perturbative expansion in an independent way, by trying to estimate 
what is the true expansion parameter for the finite-temperature theory, near the critical 
value of 0. Let us first review how the perturbative power counting argument goes in the 
standard model. The cost of adding one additional thermal loop of W bosons to an arbitrary 
Feynman diagram contributing to the free energy is 

9 2 T C r d 3 p g 2 T c Qirm w m 2 h0 ^ f m h o \ 2 

(2tt)3 J (p2 + m 2 ^ (0 C )) 2 K m w (0 c ) ~ 2m 3 w + m| " v35 GeV J ' 1 ' 

where we used mw(4>c) = g^c/^, m^o = 2A e fft> 2 , and the standard model value for critical 
ratio 4> C /T C = (/>o/T = (2m^ + m 3 z )/(2ir\v 3 ). In ref. [11] it was shown that the true loop 
expansion parameter is around 0.4 for a Higgs mass of 35 GeV, while ref. [8] finds that C /T C 
changes by 30% in going from one loop to two loops when m^o is near 35 GeV. In a two-Higgs 
doublet model, it is no longer the W boson which gives the dominant corrections to the po- 
tential, but rather the new Higgs bosons. Thus the expansion parameter will be AjT c /m 2 ((/> c ) 
instead of (17), where Aj and rrii are any of the quartic couplings and corresponding mass 
eigenvalues, respectively, of the new Higgs sector. The couplings of the mass eigenstates will 
be linear combinations of the couplings in the Lagrangian, Ai, hi, h 2 and h 3 , which are given 
by eq. (13) and the relations [4] 

hi = 2A(2yU 2 — 2m| + m 2 A o) / m 2 h o\ 

h 2 = 2A(m 2 - 2fij + 2m| - m 2 H0 ) / m 2 h o-, 

h 3 = \(2fi 2 3 -m 2 A o)/m 2 h0 . (18) 

Comparing with the standard model, it thus seems reasonable and conservative to require 
that 

e = maX ' A "^ & ' (19) 
mm(|mj(0 c )|) 

in order to have some confidence that the perturbative expansion is not out of control due 
to the infrared problem. In the last column of figure 3 we show the contours of e, using as 
inputs the thermal masses at C and T c corresponding to the third column. It is interesting 
to see that e and C /T C tend to be anticorrelated: where one is big the other is small. This 
is good news as concerns the trustworthiness of the interesting regions where C /T C > 1; it 
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0.5 0.75 1 0.5 0.75 1 0.5 0.75 1 

Figure 5: Contours of </> c /T c in the plane of m^o versus m^o, using the same three respective 
methods as in figure 3, setting the other Higgs boson masses to their tree-level values in the 
minimal supersymmetric standard model. The horizontal strip at the bottom, labeled "s," 
shows the corresponding values in the standard model. 



means we expect the value of 4> C /T C predicted by the 1-loop result to be rather stable against 
changes due to higher order contributions. 

4.1 Supersymmetry 

The model we chose to study is not quite a realistic limit of the minimal supersymmet- 
ric standard model (MSSM), even in the case that all the superpartners are much heavier 
than the Higgs bosons and therefore decouple. This is because to be consistent with the 
assumption that the symmetry breaks in the direction where the two Higgs fields have equal 
VEV's, we assumed that the top quark coupled equally to each, which is not possible in 
supersymmetry. Nevertheless one might trust the model to give a qualitative indication of 
the effects of the additional Higgs bosons on the phase transition. In the limit of tan /3 — 1 
to which we restricted ourselves in this study, the tree-level masses of the Higgs bosons are 
related according to [27] 

m 2 h o = 0; m 2 H o = m 2 A o + m|; m 2 H ± = m 2 A o + m^, (20) 

and the parameter n\ = m A0 /2, as can be seen from eq. (18) using the fact that h 3 = 
in supersymmetry. Of course the h° boson is not really massless, but it gets a mass from 
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loop corrections involving the top squarks, of order a{m 2 /mz) ln(m t /m t ). Thus we are only 
interested in rather light h° particles. For the heavier Higgs bosons on the other hand, the 
tree level relations should be a good approximation for their masses, showing that m A o is 
the only free parameter. 

In figure 5 we display the contours of <f> c /T c using the same three methods as in figure 3, 
in the plane of m A o versus m h o. For comparison, the standard model results are shown in 
the horizontal strip at the bottom of each graph. In contrast to figure 3, here we see that the 
distinction between the more accurate critical temperature and the naive one, T , is quite 
important, and furthermore the differences between the two methods of ring-improvement 
are significant. The latter observation, comparison with the standard model results, and the 
lack of variation of C /T C with m A o show that the extra Higgs bosons are generally playing 
a small role; it is the transverse gauge bosons which are most important in determining the 
strength of the phase transition. The perturbative expansion parameter corresponding to 
their contributions, eq. (17), is large for m^o > 35 GeV, hence the discrepancy between the 
two methods of resummation. Although there is a hint that very light masses for the extra 
Higgs bosons can strengthen the phase transition, the apparent instability of the perturbation 
series in this region renders any such conclusion uncertain. 

It is interesting to compare the MSSM case to that shown in figure 3. In the latter one 
sees that for a fixed value of n\ it is always possible to choose the other masses large enough 
to make C /T C > 1. But in the MSSM, there is a special relation between n\ and the new 
Higgs boson masses which prevents them from enhancing the first order nature of the phase 
transition in the minimal supersymmetric standard model, even when the bosons become 
quite heavy. This is easily understood in terms of the couplings hi in eq. (18), which are 
constrained to be of order g 2 in supersymmetry, whereas in a generic two Higgs doublet 
model they become large as the scalar masses increase. 

Also interesting is the large difference between using the naive critical temperature T 
and the more accurate one T x . The former gives a substantial overestimate of the strength 
of the phase transition. In ref. [14], regions of parameter space consistent with 4> C /T C > 1 
were found in the MSSM with tan (3 = 1, however using T as the critical temperature. Our 
results indicate that one should use a more precise definition of T c before concluding that 
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the Higgs bosons can have much effect on the phase transition. 

5 Conclusions 

We have investigated the electroweak phase transition in a class of models containing two 
Higgs doublets, with a global symmetry that insures electroweak symmetry breaking along 
the direction $1 = $ 2 . The results are encouraging from the standpoint of electroweak 
baryogenesis, which requires that <j) c /T c > 1 to avoid the washout by sphalerons of the 
baryons that might be produced. Even for rather heavy h° Higgs particles up to 300 GeV, 
it appears possible to choose the others sufficiently heavy to make a strongly first order 
transition, except perhaps when the parameter /x| from the mass term /i|$|$ 2 becomes too 
large. In supersymmetric models however, /z§ is automatically tuned in such a way as to 
minimize the impact of the extra bosons on the strength of the phase transition, no matter 
how heavy they are. 

We tried to make our study as reliable as possible by using certain improvements: in- 
clusion of Goldstone boson contributions, going beyond the small M/T expansion, and esti- 
mating the size of higher order perturbative effects by comparing different methods of ring 
diagram resummation and by computing the expansion parameter of the finite temperature 
theory. The latter checks give us confidence that our conclusions are quantitatively correct. 
However it is always possible that perturbation theory fails despite such safeguards. Recently 
it has become possible to start investigating the electroweak phase transition in theories with 
new physics using nonperturbative lattice results, by mapping the theory of interest onto 
an effective lagrangian of the same light degrees of freedom as in the standard model [28], 
which has already been analyzed on the lattice. Steps in this direction were recently taken 
for two Higgs doublet models by reference [29]. It will be interesting to compare the nonper- 
turbative approach with that pursued here, to verify more conclusively that there really are 
conditions when perturbation theory, which is easier to apply and more amenable to analytic 
expressions, can be trusted to give quantitatively accurate results for the finite temperature 
effective potential. 
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